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Abstract 

We use the method of finite temperature field theory to examine 
the propagation of mixed neutrinos through dense media, putting the 
emphasis in those situations in which the neutrinos themselves are 
in the background. The evolution equation for the flavor amplitudes 
is deduced, and the expressions for the corresponding hamiltonian 
matrix are given explicitly. We find that, in order to include the non- 
linear effects due to the I'-i^ interactions, the neutrino propagator that 
must be used in the calculation of the neutrino self-energy diagrams 
that contain neutrinos in the internal lines is the propagator for the 
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neutrino modes in the medium instead of the thermal free-field prop- 
agator. We also show how the absorptive contributions are included 
in terms of a non-hermitian part of the hamiltonian matrix, which 
we indicate how it is calculated. Our treatment provides a consistent 
generalization of a method that has been successfully applied to the 
study of neutrino oscillations in matter. 

1 Introduction and Summary 

Finite Temperature Field Theory (FTFT) provides a very useful framework 
to study the behavior of elementary particles, including neutrinos, in the 
presence of a thermal background|T[]. A quantity of fundamental interest 
in this formalism is the self-energy of a particle, from which the dispersion 
relation of the modes that propagate in the medium are obtained. In the 
case of neutrinos, the density and temperature dependent contributions to 
the self-energy can be expressed as an effective potential for each neutrino 
flavor that, added to the vacuum kinetic energy, yield the Hamiltonian ma- 
trix that governs the neutrino oscillations in matter p|, ^. In normal matter 
(electrons and nucleons) the universality of the neutral-current interactions 
of the neutrinos implies that their contribution to the effective potential is the 
same for all flavors and hence irrelevant as far as oscillations are concerned. 
However, in environments like the early universe or the core of a supernova, 
where the neutrinos represent an appreciable fraction of the total density, 
the contributions to the potential energy arising from the v-v scattering are 
not in general proportional to the unit matrix and should be included in the 
analysis of the resonant flavor transformations. 

The calculation of the contribution from the background neutrinos to 
the effective potential presents a new challenge because now the background 
depends itself on the flavor amplitudes, whose time evolution is in turn deter- 
mined by the oscillation mechanism. Thus the problem becomes a non-linear 
one that must be solved in a self- consistent manner. Several studies of the as- 
trophysical and cosmological implications of neutrino oscillations under these 
conditions have been carried out during the last years[^, ^, These works 
left out the off-diagonal contributions of the v — v scattering to the flavor 
energy matrix. The existence of these extra terms was first recognized by 
Pantaleone[0] and later on have been included in certain numerical analysis 
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of neutrinos oscillations 

Most of the extensive studies of neutrino oscillations in matter and the 
MSW solution to the solar neutrino problem|T^ are based on the Wolfenstein 
equation for the flavor amplitudes of a single-particle wave function. On the 
other hand, according to the prevailing attitude in the literature on the sub- 
ject, a consistent description of neutrino oscillations that includes the effects 
of their mutual interactions and the noncoherent interactions with the back- 



ground is possible only within the density matrix formalism||g, |TT], |T2|, |13 
The works cited in Ref. |jl2|, [T^ give a general treatment of relativistic mixed 
neutrinos along this line which include the neutrino self-interactions, and 
have been recently applied to examine the matter-enhanced flavor transfor- 
mations in supernovae[|I^, |T^]. 

While it is widely accepted that FTFT is an efficient method to compute 
the background effects that are incorporated as a potential energy in the 
Wolfenstein equation, a description based on this approach when neutrinos 
compose the background is lacking or incomplete. The present work fills this 
gap. By working within the formalism of the FTFT, we derive the evolution 
equation for the fiavor amplitudes of a neutrino that propagates in a medium 
containing several species of massive neutrinos. Explicit formulas for the 
energy matrix are given in a familiar form that make them amenable for 
application to concrete situations. In order to bring out the salient features 
without introducing unnecessary complications, we work in detail the case of 
two Dirac neutrinos in the Standard Model, but our approach makes evident 
the path to treat more general cases, such as including additional families 
and Majorana neutrinos. 

One of our objectives is to show that FTFT, which is simpler to use than 
other approaches, can be consistently applied to examine the type of problem 
we have referred to. The procedure is to calculate the self-energy matrix for 
mixed neutrinos, which leads to the effective Dirac equation obeyed by the 
modes that propagate in the medium. We obtain the Hamiltonian matrix 
for the flavor evolution equation directly from the condition that determines 
the energy-momentum relation for the neutrino modes. A particular new 
feature that arises when neutrinos appear in the background is that the 
neutrino densities are modulated by the oscillation mechanism itself, but on 
the other hand the standard FTFT propagators do not take this into account. 
Therefore an essential part of the problem is to flnd the appropriate neutrino 
propagator to be used in the self-energy diagrams that contain neutrinos 
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in the internal lines. This aspect of the problem is considered in Section ^ 
where we show that this is not the standard free-particle thermal propagator, 
but rather the one for the effective field constructed out of the one-particle 
wave functions of the propagating modes. The use of this propagator is a key 
ingredient in our treatment. We also show how to include the damping effects 
by means of a non-hermitian part in the Hamiltonian0], which is determined 
in our scheme from the absorptive part of the self-energy. In Section ^ we 
prove that the non-hermitian part of the Hamiltonian can be expressed in 
terms of the total rates for the emission and absorption processes of the mixed 
neutrinos by the background. There we also discuss some subtleties that arise 
if these rates are determined in the traditional way, in terms of probability 
amplitudes, when neutrinos compose the background and we explain how 
they are overcome in our treatment. 



2 Normal Matter 

The dispersion relations and wave functions of the neutrino modes that prop- 
agate through a medium are determined from the linear part of the effective 
field equation. In momentum space it takes the form 

(^-m-Se//)^ = 0, (1) 

where fc^ is the momentum vector and m is the diagonal matrix whose ele- 
ments are the neutrino masses in vacuum. In the basis of the mass eigenstates 
z/j, the background contributions to the neutrino self-energy Se// is a non- 
diagonal matrix which can be written as[|l^ 

{^eff)ij = {dLijf^ + hujii + CLij)L + {anijf: + hmjii + CRij)R , (2) 

where is the velocity four-vector of the medium and L,R = (1 =p 7^)/2. 
In the rest frame of the medium 

(i,o) , 

and the components of are given by 

P = (w, k) . 
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In general, the elements of the matrices aL,R, bL,R and cl,r in Eq. (|^) 
depend of u and but in a isotropic medium they are functions only of the 
Lorentz scalars 

u = k ■ u , 

K = VUJ^ - fc2 , (3) 

with K = |k| . In what follows, the arguments indicating the dependence on 
these variables will be generally omitted, but they will be explicitly included 
when needed for clarity. 

The background is taken to consist of nucleons, electrons and neutrinos, 
and their respective antiparticles. We assume that the temperature is low 
enough so that the contributions of the heavier charged leptons can be ne- 
glected. The relevant piece of the weak interaction Lagrangian is 
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■ i f=e,n,p 



^W^Y.Ue^eLl,'yL^, (4) 

where the index i runs over all the neutrino species. The un are the left- 
handed components of the neutrino fields with a definite mass, and U is 
the unitary matrix that relates them with the neutrino flavor fields z/^q, = 
J2iUai^Li (tt = r). The coefficients Xf and Yf are 

Xe = -^ + 2sm'^ew, 

Ye = I, (5) 



for the electron and 



2 
1 



Xp = - - 2 sin 6*1^ , 
_ _ _ 1 

Yn = = —Yp = 2 ' (^) 



for the nucleons. In addition, 



mz cos 9w = TTiw , 

= V2Gf. (7) 



9 
Aw? 
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The calculation of the self-energy proceeds as in the vacuum, but with 
the free propagators for the internal lines replaced by their thermal general- 
izations. Since we are assuming that there are no W and Z bosons in the 
medium, their propagators are the same as in the vacuum. 

Let us first discuss the simplest case of a homogeneous background that 
does not contain neutrinos. For Dirac neutrinos propagating in such a back- 
ground, the result of the one-loop calculation of the self-energy, up to terms 



of order g /m^, is|T7 



aL,Rij — CL,Rij — , 

bRij = , 

biij = bij = V2GF[U*iUej{ne-ne) +6ijQz] , (8) 



where 



Qz= E Xfirif-nj), (9) 

f=e,n,p 

is the average Z-charge of the medium. In these formulas, Uf and rij denote 
the total number densities of the particles (electrons, neutrons and protons) 
and the antiparticles, respectively. 

In terms of the left- and right-handed components of ip Eq- (0) becomes 

4Li^L-m^R =0, 

jf^^iljR-miPL =0, (10) 

where 

-A-^L = — bUfj, , 

A^R = k^. (11) 

In writing these equations we have used the particular results given in Eq. (|^). 
In a more general case A^l r = (1 — a^ Rjk^ — b^RU^ and the matrix m in 
Eq. (p!OD has to be replaced by m+CR^L in the first and in the second equation, 
respectively. 

We use the Weyl representation of the gamma matrices and put 

^R = 

= ( : ) . (12) 
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The equations to be solved then become 



(13) 



Since A^ and Aji are both proportional to k, these equations can have non- 
trivial solutions only if ^ and t] are proportional to the same spinor (px with 
definite helicity, defined by 



with A = ±1. Therefore, we write 



(14) 



V 



(15) 



where y and x are spinors in the basis of the z/j. For a given helicity A, the 
equations for y and x that follow from Eq. (|1^) are 



A\ X — my 
A\j^ y — mx 



0: 

0, 



(16) 



where 



4^^ 



A 



(A) 
R 



uj + \k — b . 

uj + Xk, . 



(17) 



Using the second equation in (|T6D to eliminate y from the first one, the 
equations for y and x become 



A'-^^'^m-^A'-L^ - m 



X 



0, 



A 



(-A) 
R 



mx 



which have solutions only if 

Det (^A'-^^^m-^A 
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m 



0. 



(18) 



(19) 



The roots of this equation determine the dispersion relations of the prop- 
agating modes. The complete characterization of the different solutions re- 
quires, in addition, the determination of the wave function for each mode. 
For these purposes it is instructive to notice that for one (unmixed) Dirac 
neutrino, Eq. ([T9|) reduces to 



K 



h){uj 



{uJ + K — b){uj — K,) 



m 



m 



(20) 
(21) 



where the expression for h is obtained from Eq. (|]) by making obvious 
simplifications. For example, for an electron neutrino in normal matter, 
h = \/2GF{ne + HfXfUf). Each one of Eqs. ( |20| ) and (pl|) yields two solu- 
tions: one with a positive real part of uo and another with a negative real part. 
The four solutions correspond to the positive and negative helicity states of 
the neutrino and the antineutrino. In particular, the dispersion relation for 
the left-handed neutrino corresponds to the positive solution of Eq. 



1/2 



+ 



(22) 



For future purposes it is useful to notice here that, for relativistic neutrinos, 
the effect on the dispersion relations of a non-zero value of the coefficients aj^ 
and a^j is negligible. To see this notice, for example, that if a^, 7^ Eq. (^) 
is modified to 



[(1 - aj}j{ijj - k) - 6](a; + k) 



m 



(23) 

For cij ^ K, the corrections to the solution given in Eq. ( |22D are either of 
higher order in or of order vp?g^ j k and we ignore them. 

Based on these remarks, let us now return to Eq. (|1^) and look for solu- 
tions with A = —1. As already mentioned, for clarity we consider the case 
of mixing between only two Dirac neutrinos. The condition that determines 
the energy-momentum relation is given explicitly by 



UJ 



K 



2 2 

UJ — K, 



where 



bp cos^ 9 



- (cj + K)b22 - ml 

{uJ + K)^bi2b21 = , 



V2GpQ2 



(24) 
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622 = heSm^e + V2GFQz 
bi2 = ^21 = be sin 9 cos 9 . 

Here 9 is the vacuum mixing angle and 

be = V2GF(ne -Ue) . 



(25) 



(26) 



Eq. (|24D has in general two positive and two negative solutions, the latter of 
which correspond to the antineutrinos. We seek approximate solutions to this 
equation in the relativistic limit, which is the adequate regime for neutrino 
oscillations. We do it first for neutrinos (positive solutions), and afterwards 
apply the same considerations to the antineutrinos. Dividing Eq. (^^ by 
{uj + k)"^ and using the approximation 



UJ + K 



which is valid for mi2, hj, Eq. 



2k ' 

reduces to 



(27) 



UJ — K 



mi 
2k 



UJ — K 



J22 



7122 

2k 



12(^21 



with solutions 



m? + ml bii + b 
uJi,2 = n-\ : h 



4k 



22 ^ 

2 ^2 



bu - &22 + 



m^ 



2k 



+ 46l2&21 



(28) 



1/2 



(29) 

With each dispersion relation uji^2 there is associated a two-component vector 
X = ei^2, which is obtained by solving 



A^^m ^ - m 



The corresponding y-vector is given by 

m 



'^='^1,2 



ei2 = 0. 



yi,2 



A 



(+) 

R 



ei,2 



(30) 



(31) 



and together with Eq. (|30|) , it determine the wave function of the propagating 
mode. Clearly, for relativistic neutrinos, yi^2 ~ (^/2K)ei,2 • 
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Using the same approximation of (|27|), Eq. (|30| ) can be recast in the form 

Hci = uJiCi , (32) 

where 



2 

m 



H = K + — + b{K, k) . (33) 
2k 



For the two-flavor mixing case we are considering, 



(34) 



with 

U-i ""flV (35) 

\ — sm u cos y / 

The corresponding Hamiltonian in the flavor basis is 

Hfiav = UHU^ , (36) 

and after subtracting an irrelevant term proportional to the unit matrix, it 
can be written in the form 

rr _ 1 / -Aocos2^ + 6e Aq siu 20 \ , . 

^'■^^ 2 1^ Aosin20 Aocos20-6ey' 

where 

A„ = ^^. (38) 

This is the energy matrix traditionally used in the studies of neutrino oscil- 
lations inside the Sun. Its eigenvalues are simply =f|(^2 — ^^i), with 



be 



002- UJl = y ( Ao COS 26 - K) + Al sin' 26 . (39) 
The explicit solutions of Eq. ( ^0|) in the flavor basis are readily found to 

^ / cos6'm \ 
[-sm6^)' 
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where 6m is the mixing angle in the medium given by 
• oQ Aosin2^ 

sm26'm = , . (41) 

'(Ao cos sin 26 



Written as in Eq. (0), the dispersion relations 001^2 are simply the eigen- 
values of the Hamiltonian H with the being the corresponding eigenvec- 
tors. The physical interpretation that emerges is the following. The Dirac 
wave function for a relativistic left-handed neutrino with momentum k that 
propagates through a medium with a definite dispersion relation uji is 

= e'"^-"^ ^ e-'^'^'ei , (42) 

where we neglect ipR, which is of order m/2K as compared to the left-handed 
component. If the neutrino is not initially in a state corresponding to one of 
the eigenmodes, then the wave function is 

i'L = e^'-' ( ^_ ) X{t) , (43) 
where xit) is the general solution of 

4^//X-. (44) 

that satisfies the specified boundary condition for x(0), and can be written 

as 

Xit)= E(elx(0))e--*e,. (45) 

i=l,2 

Since neutrinos are produced in states of definite flavor, in neutrino oscilla- 
tions problems the initial conditions are commonly specified by giving the 
flavor components of the initial state. Thus, if we denote 

X(0) = ( ; ) . (46) 

in the flavor basis, then at a later time t 

xit) = x^'\t)ae + x^'^Ht)a,, (47) 
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where x'^'^'^H't) the solutions to Eq. ( PD with the initial conditions 



1 







1 



(48) 



Physically, Og,^ is the amplitude for finding the neutrino initially in the state 
corresponding to jt'c/i)- If we write 



x{t) 



Xe{t) 



(49) 



then 



Xa (t) = tteX'-a^ i^) + ^f^Xa^ (^) , (« = 6, /i) 

is the amplitude for finding the neutrino in the corresponding fiavor state at 
subsequent times, with a similar interpretation for the components of x{t) in 
other basis. 

Associated with the wave functions there is a quantum field 







ix^'\t)a,{K) + x^'^\t)a,iK)) 



(50) 



where x^^'f^'i are the functions defined above, while the ae^^{K) are now in- 
terpreted as the annihilation operators of left-handed neutrinos of definite 
fiavor with momentum k. 

Up to now we have restricted ourselves to the positive frequency solutions 
of Eq. ([T9| ) with A = — 1. In order to include the right-handed antineutrinos 
we need to consider the negative solutions with A = +1. Writing the negative 
solutions as a; = — oJ, then uji^2 are given by the same formulas as in Eq. (p^) 
but with the sign of the coefficients hij reversed: 



_ m? + fell + 622 1 

^^'^ = '^+^;^ 



J22 



bii + 



mi — 
2k 



+ 4612621 



1/2 



(51) 



The equations corresponding to Eq. (|30D for the two-component vectors x 
ei 2 are 



- ei,2 







(52) 
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Similarly to Eq. (p^, they can be recast in the form 

H_ei = uJiEi , (53) 

where 



2 

m 



H_ = K + — - b{-K, k) . (54) 



In particular, if_ differs from Eq. ( plf ) only in the sign of the Qz and be 
terms. In analogy with the positive frequency field \ we built a negative 
frequency field 



where the aeji{K) represent the annihilation operators of (right-handed) an- 
tineutrinos of definite flavor, with momentum k. The x^^'^'^ solutions 
of 

4=7rx, (56) 

that satisfy 



1 




X^^HO) = ( ? J , (57) 



with 



2 

m 



K + - bl{-K,-K) 

2k, 



The corresponding matrix in the flavor basis is 

Hflav = U*HU^ 

_ 1 / - Ao cos 26 - be Ao sin 26 
~ 2 V ^0 sin 26 Aq cos 26 + b^ 



(59) 
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where we have again discarded a muhiple of the identity. Two comments 
are in order. Although we are considering a homogeneous medium, in some 
of the previous expressions we have retained the dependence in k to make 
more obvious their generahzation to the anisotropic situation. In addition, 
for the case of mixing between two Dirac neutrinos, the matrix U is real (see 
Eq. (|35|)), but in the case of Majorana neutrinos three Dirac neutrinos, this 
is not not necessarily true; for that reason, in Eq. (^) we have not replaced 
U* by U. 

Neglecting the right-handed components then, the complete effective neu- 
trino field is 

^iJ{x) = ¥^\x) + ¥p{x). (60) 

The wave functions used to define "^{x) must include appropriate normaliza- 
tion factors Nf^. They can be calculated in terms of the self-energy [p!^] and, 
in general, differ from the normalizations in the vacuum. With the coeffi- 
cients as given in Eq. (H), N^^ = 1 and therefore, if we take the creation and 
annihilation operators normalized such that 

{al,{K'),aa{K)} = {27Tf6^^\K' -K)6a'a, (61) 

then the spin wave functions 0a have to be chosen to satisfy 

0V = 1 • (62) 

So far we have assumed that the background is homogeneous. Now sup- 
pose that the matter density is not constant. Then along the neutrino path 
the densities of the different particles in the background become functions 
of the time t. The crucial step that we take is to postulate that, in such 
situation, Eq. (^31) and (|^) remain valid and determine the spinor wave 
functions of mixed left-handed neutrinos with momentum k. Consequently, 
Eq. (|50|) is assumed to give the associated effective field also for an inho- 
mogeneous medium. Of course, in that case Eq. (^) no longer holds, and 
an appropriate solution of Eq. (|4^) must be found. Identical considerations 
applies to the (right-handed) antineutrinos. Our Eq. (^if ), with the (time- 
dependent) Hamiltonian expressed in terms of be{t) = \/2GFne{t), is pre- 
cisely the Wolfenstein equation that determines the flavor neutrino evolution 
through ordinary matter 0. 
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3 Neutrino background 

As illustrated in the previous section, our formalism reproduces the standard 
formulas that have been utilized in the numerous studies of the MSW effect 
when it is applied to a medium made of electrons and nucleons. In what 
follows we show how to extend this field-theoretic approach to take into ac- 
count the new features that appear when several species of neutrinos are also 
part of the background. In this situation, the neutrino self-energy receives 
additional contributions from the diagrams depicted in Fig. 1, which involve 
neutrinos propagating in the internal lines of the loop. These contributions 
are given by 

(63) 

X (-l)Tl'/ ^ fei^"' (*')«) T-i. (64) 

where the trace is taken over the Dirac indices. S^p\k') is the thermal neu- 
trino propagator defined by 

' I = (^^^(^)^.(?/)) . (65) 

where T stands for the time-ordered product of the fields. The angle brackets 
in Eq. ( |65D indicate a thermal average, defined by 



for any operator (9, where 



and 



{O) = ^Y.^^\pO\^) , 

p — q-I^'H+Y.a^aQa 
i 



(66) 

(67) 
(68) 



In Eq. (^) (3 is the inverse temperature, Ti is the Hamiltonian of the system, 
the Qa are the (conserved) charges that commute with and a a are the 
chemical potentials that characterize the composition of the background. 

In the ordinary perturbative applications of FTFT, the right-hand side 
of Eq. (|65|) is evaluated in terms of the free fields. Then, a straightforward 
calculation of the thermal averages of the products of creation and anni- 
hilation operators yields the canonical formulas for Sp{k) in terms of the 
particle number densities. In particular, the propagators for the electron 
and the nucleons calculated that way, were used to compute the contribu- 
tions to the neutrino self-energy given in Eq. (|^). However, for the problem 
at hand we cannot use the thermal free neutrino propagator to calculate the 
additional contributions to the self-energy introduced by the neutrino mu- 
tual interactions. The reason is that the neutrino densities are modulated by 
the oscillation amplitudes, which evolve according to the equations we have 
obtained (Eqs. ( PI) and (p3D), and the thermal free neutrino propagator does 
not incorporate this effect. 

Our proposal in this work is to use, instead of the thermal propagator 
of the neutrino free-field, the propagator that corresponds to the effective 
neutrino field defined in Eq. (|60D. The calculation of this propagator proceeds 
as in the free field case, simply by substituting the expansions in Eqs. (^) 
and (|56| ) into (pS)). However, as will be seen shortly, for our purpose here, 
it is not necessary to compute it explicitly. The result of a straightforward 
calculation, whose details we give below, are 



v^G^(7^L)(j>)+jf ),,, (69) 
5,,^Gf{^^L) J l^n; [U - fv. + (e - /x)] , (70) 



where 



K 



(27r)3'> 



X^'h^'^^U. + (e - /i) 



(27r)3 ^ 

The vector n'^ has the components 



ix^V^^Th. + (e - /i) 



(71) 



n 



(72) 
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in the rest frame of the medium, and f^^ {f-p^ ) are the momentum distri- 
butions, normahzed such that 

J (27r)3^"^''^' 



are the number densities of the flavor neutrinos (antineutrinos) at the initial 
time. 

Let us first derive the result for T.^'"^\ Using the local approximation for 
the Z-boson propagator, 

AfJ^^, (74) 

and the Fierz-like identity 

7"LA7„L = -(TrA7,L)7"L, 
which is valid for any 4x4 matrix A, we have 

(A:) = -^ i^L) J ^Jr{SAk%rL) . (75) 



As in Eq. (^^, the trace is with respect the Dirac indices. While we can 
calculate the propagator using the definition in Eq. (|65|) and then substitute 
the result in the above formula, it is more expedient to observe that 

- ^ / ^Tr(5^(A;').,7^^) = (iV^,(x)7,^.(x)) , (76) 

which results immediately from Eq. (|65D , by taking the limit x — > In the 
last equation N stands for the normal-ordered product. The quantity in the 
right-hand side is easily computed by inserting the plane wave expansions 
given in Eqs. (pDj ) and (|^). The statistical averages are evaluated by means 
of the formulas 

(a^,(K')aa(K)) = i2nf5^'\K' ~ K)5^'afv^ (77) 
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where a = e, jj,. In this manner we obtain 



(iVvl/^-(x)7,vl/,(x)) = (j>)+jf)- 



(78) 



and using Eqs. (0) and (|7|) in (|7|) we finally arrive at Eq. (|9D. For ^'^""^ 
by a similar procedure it follows that 



Since the hamiltonian matrix is hermitian, 
and 

^ + jP)u = J [/- - + (e - ^)] 



which gives Eq. (|70|). 

T.^'^"^ can be decomposed as in Eq. 



with 



V2G 



X 



(27r)3 V / 



K 



K ■ K 
1 U 

K 



x^'h^'^^U - {x^h^^^Th. + (e /i)] 



(79) 

(80) 
(81) 



(82) 



(83) 



As remarked below Eq. (^2|), the coefficient a'^'^'^^ has a negligible effect on the 
dispersion relation and therefore we do not consider it further. From Eq. (^) 
and the results of Ref. |18], it can also be verified that the wave function nor- 
malization to which we alluded before remains equal to one. On the contrary, 
}){'^'^) produces a new significant contribution that modifies the hamiltonian 
H given in Section 0. Since, for relativistic neutrinos (antineutrinos), u ^ k 
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(to ^ —k) the above formula for h^'^"'' translates into the following additional 
term to the matrices exhibited in Eqs. ( ^If ) and (|58|) : 

-jji'^-) = = (84) 

The additional terms produced by S''^'^'^^ can be obtained in similar fashion; 
they are proportional to the identity matrix I and are given by 

^,[uu) ^ _;^'(-) ^ ^CpQ^^r'^I , (85) 

where 

^^''^ = 10)''^^'^''''^ ~ ^""^ ■ ^^^^ 

To summarize, the propagation of mixed neutrinos through a background 
that includes neutrinos is described by the effective neutrino field given 
in Eq. (|60|). The time evolution of the amplitudes x^"\'t) ^^id x^'^H'^) that 
appear in the Fourier expansion of the field is governed by 

I ^ = (^]^{rnatter) _|_ j^{vu) _|_ JJl{uu)■^^{e,^J.) ^ 

dt 

where and f/'*^""^**'^'') stand for the normal-matter contributions de- 

termined in the previous section (Eqs. (|3^ ) and (H)), with the label matter 
added to single them out. Discarding a term proportional to the identity 
matrix, in the flavor basis the above equations can be put in the form 

■ dX^ = 1 / -AoCOs20 + 6e + 2/lee Aq siu 2^ + \ (,.^) 

^ dt 2\ Aosin20 + 2V Aq cos 2^ - 6e + 2 V / ^ 

.dx^ = 1 f -^0 cos 26 -K- 2/iee Ao sin 26 - 2K^ \ (^,^) 
' dt 2 I Ao sin20 - 2V Aq cos 2^ + K ~ 2h^^ ) ^ ' 



(88) 



where h^n' are the elements of 



2 
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and are explicitly given by 



-h 



V2Gt 



(27r)3 

Xe 



K ■ K] 



a=e,fj. 



''e/j. 



h 



(27r)3^ 



I Xl 



K ■ K] 



fu. 



\xll I ) fua 



fu 



(a) (a)*\ r_ 
Ae A/t I J Ua 



(89) 



(90) 



a=e,fi 



In the present case, in which we are neglecting absorptive effects, we can use 



the probability conservation conditions 
and get 

d^K' 



1 and Ix^^^P = 1 in Eq. (||) 



he 



-h 



MM 



V2Gj 



'1 



K ■ K] 



a=e,fi 



x':'\'- 



1 



fu 



(91) 



The set of (non-linear) equations in (|88|) (or (q^)); subject to the initial 
conditions given in Eq. ( ^8]) and (|57D, are the ones that must be solved 
in the application to an actual physical problem [|l9l. They represent the 
generalization of the Wolfenstein equation to the physical situation under 
consideration. In our notation, hf,^^ and h^^e account for the nondiagonal 
contributions to the potential energy in the flavor basis, that arise from the 
z/z/ interactions 0. 

In deriving the previous results no assumption has been made on the 
characteristics of the neutrino background. If the medium is isotropic, then 
the quantities between brackets in, for example, Eq. ( P5| ) can not depend on 
the angular integration variables. In such the integral involving the 

factor k' ■ k vanishes and the coefficients in the self-energy decomposition 
reduce to 



0, 

V2Gf 



d^ 



K' 



{2nf 



X^^^^'^^fu. - {x^^h^^^Tfv. + (e - /i) 



(92) 
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with similar simplifications in Eqs. ( |89| ) and (pOl). 

Usually, the quantity of interest in applications is the density of electron 
neutrinos at a given time. In our language it is simply 



Ml 



(93) 



which, upon substituting the expansion given in Eq. (|50|) and evaluating the 
statistical averages by means of Eq. (^), yields 



K 



(2vr)= 



\X. 



(e)|2 



(94) 



The density of neutrinos of other types can be computed in identical manner 
in terms of the initial distributions and the flavor amplitudes. For example. 



(27r)= 



\X 



(e)|2 



|y(M)|2 



(95) 



with analogous expressions for antineutrinos. At t = 0, these formulas for 
and n,^^ reduce to their appropriate initial values defined in Eq. ([731). 
At each time t, the Hamiltonian matrix in Eq. (|88|) can be diagonalized 

by the unitary transformation Um{t) 











cos6'r, 



sin 6r 



sin 9rr, cos 9r, 



(96) 



with the time-dependent mixing angles 9m and defined by 

sin 2»„ = + . (97) 

',9 -be- 2heef + (Ao sin 29 + 2hef,Y 

arg[Aosin2^ + 2/ie;,]. (98) 

Under the present conditions, the difference between the (instantaneous) 
energy eigenvalues is given by 



U!2 — UJl = y ( Ao COS 29 — be — 2he 



(Ao sm29 + 2he^)^. 



(99) 
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As is evident from Eq. (^) the neutrino mixing angle in matter is modi- 
fied by their mutual interactions, and the resonance condition (sin26'm = 1) 
becomes 

Ao cos 20 = 6e + 2/lee • (100) 

Moreover, the off-diagonal elements of the Hamiltonian are in general com- 
plex, and the phase will be present in contrast with the ordinary-matter 
situation. Similarly, a unitary matrix Um{t) can be introduced for antineu- 
trinos; the respective formulas for the angles 9m and are obtained by 
replacing be, hee, and hefj, by their negatives in Eqs. (|97|) and (p^). 



4 Absorptive effects 

In the previous sections we have assumed that the corrections of order g"^ /m^ 
to the real part of the self-energy are not important and in the same spirit 
we have discarded the imaginary part of the self-energy. This is normally 
justified since the corrections of order which are proportional to the 

particle-antiparticle asymmetries, give the dominant contributions. However, 
this is not true for a CP-symmetric plasma like the early Universe, and in 
such situations the real and the imaginary part of the neutrino potential 
become comparable |2^. Our purpose now is to complement the formalism 
by examining the effects induced when the imaginary part of Se// is not 
negligible. 

The damping effects, which are due to the incoherent interactions of the 
neutrinos with the particles in the background, are described through an anti- 
hermitian part of the Hamiltonian that governs the evolution of the flavor 
amplitudes. In our treatment, the damping terms in the Hamiltonian are 
determined from the imaginary part of Sg// in a way that we now explain. 

We have learned that the complete Hamiltonian for neutrinos is 

2 

TfL 

H = K + — + 6(k,k), (101) 

where b{uj, k) stands for the coefficient of in the decomposition of Sg/j in 
Eq. (H) . In terms of the vector defined in Eq. ([72|) , we have 



h{K,K) = ^(TrLjiSe//) 
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= ^4Se//(«,/?)H°, (102) 

where, in the second equahty, we have introduced the massless vacuum Dirac 
spinors normahzed such that 

ulul = KLfi. (103) 
If Se// develops an absorptive part, so that 

7°St^^7VSe//, (104) 
the matrix b develops a non-hermitian term. Thus, if we write 

Ee// = + iEj (105) 
b = br + ibi (106) 

where br,i are hermitian matrices and 

1 

s< = i 

2i 

are the dispersive and absorptive part of the self-energy, then 



■i^eff-l'^lffl') (107) 



b, = i(TrL^EOL=, 

= ^II°E,(«,«X. (108) 

This formula allows us to calculate the absorptive terms in the evolution 
equation for the amphtudes x^"\ from an exphcit calculation of the absorp- 
tive part of the self-energy. However, using a very general argument, we will 
now show that the matrix bi is simply related to the total rates for decay and 
inverse decays of the neutrinos. 

The starting point of our discussion is the formal definition of the elements 
of the self-energy matrix in the real-time formulation of FTFT 

i^2i{z -y)ai3,ij = -{ria,i{z)r]pj{y)) , 

i^uiz -y)ai3,ij = {r]f^j{y)r]a,i{z)) , 

-J^iiiz-y) = j:,,{z-y)e{z'-y') + J:,2{z-y)e{y'-z'), 

-T.22{z-y) = E2i(^-|/)^(/-^°) + Ei2(^-Z/)^(^°-/),(109) 
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where a, (3 are Dirac indices and i,j are indices labeling the neutrino mass 
eigenfields. rji and rj^ are the neutrino source fields, and in terms of them the 
interaction Lagrangian is 

^int = ^lV + V^L ■ (110) 
In terms of the elements Hab the physical self-energy is given by [pi] 



Se// = Sn+Sl2, (111) 

which from Eq. ( p.09| ) is immediately recognized to be the retarded self-energy. 
Decomposing En in analogy to Eqs. ( |105| ) and ( p.07] ), Eq. ( |111| ) is equivalent 
to 

Ej.(fc) = Eiir , 

Si(fc) = Sin-^Si2(A;). (112) 

On the other hand, using the integral representation of the step function, it 
follows from Eq. ([109D 

^m{k)='-[^2i{k) + ^u{k)] , (113) 



which implies that Ej, determined from Eq. (112), can be equivalently com- 
puted using 

S,(A;) = ^[S2i(A:)-Si2(A:)] . (114) 

Then from Eq. (|T08D 

&, = £-w° [S2i(fc)-Si2(A;)]w°, (115) 

which is a useful formula because the matrix elements of S21 and S12 are 
related to the rates for emission and absorption of the neutrino. To see this, 
we insert a complete set of states between the field operators 77 and in 
Eq. (|T09|), giving 

n,m 

(nh(0)|m)(m|r7(0)|n)Z„, (116) 
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where Z is defined in Eq. (p^), 

Zn = {n\p\n) 

= e-^'?"-"+"", (117) 

and an is the eigenvalue of the operator a = J2a c^aQa corresponding to the 
state |n); i.e., a\n) = an\n). From Eq. ( |116|) we immediately obtain 

iulJ:2i{n, k)uI = -2kTd (118) 

where 

= ^^E(^I^L^(0)l"^)("^l^(0)<|ri)(27r)V^)(A; + g„-gjZ„. (119) 

In particular, the diagonal elements {T£,)ii are the total rates for the processes 
n+Ui —>■ m, averaged over all possible initial states. The off-diagonal elements 
have a similar structure but they involve the product of the amplitudes for 
two different neutrino species. In similar fashion, 

ilfL^l2{K,K)ul = 2KTi (120) 

where F/ is the matrix 

= ^^Y.(M<Vmm){m\r]{0y»{27r)'6'-^\k + g„ - qjZ^ , (121) 

^ ^'^ n,m 

whose diagonal elements are the total rates for the inverse processes m 
n + Ui, averaged over the initial states. Then using the results of Eqs. ( |11^ ) 
and ( |120| ) in Eq. ( |1 1 5| ) we finally arrive at 



b^ = ~, (122) 

with 

F = Fz5 + F7. (123) 

In situations of exact equilibrium, the share the same chemical poten- 
tial a^. In that case, since rj has the same quantum numbers as the neutrino 
field, then 

{n\u{p)ri{0)\m) ^ 
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only for those states such that am = an + oiy . Thus, in Eq. ( |121| ) we can 
replace 

Z^ = e-^«+"''Z„, (124) 

which yields 

r = (l + e-^)rB, (125) 

with X = (3n — a^. However in situations that do not correspond to exact 
equilibrium the formula in Eq. (|123|) is the appropriate one. 

Then in summary, if we denote by the energy matrix displayed in 
Eq. (1^ 

H,. = + ifC^'-) + H'^'^"^ , (126) 

then the complete Hamiltonian describing the attenuation effects induced by 
the incoherent neutrino interactions is 

H = Hr-t^. (127) 

We remark that, since the total probability is not conserved any more, the 
replacements 

I (a,a)|2 ^ 1 _ I (a,a)|2 
lA/i I ^ lAe I 

made in Eq. ( pi] ) are not valid, and the correct formulas for h^.^. and /i^^ are 
those given in Eq. (|89D in this case. 

For the antineutrinos, the expression equivalent to Eq. (|101D is 

2 

ryi 

H = K + — -h\-K,-K), (128) 

which determines the evolution of the amplitudes x*^""*- Since 
h*{-K,-K) = ^ (TrL^Se//(-K, -k)) 

= j-if^E,ff(-n,-K)ul, (129) 
it follows that the complete Hamiltonian becomes 

H = H,-t^, (130) 



26 



where 



r = Td + Ti. (131) 



Tdj are given by the same formulas for F/j / in Eqs. ( |119| ) and ( |121| ), but with 



the substitution k — >■ —k in the delta function. By crossing then, it follows 
that the diagonal elements of To and Tj are the rates for the processes 
n — > m + I7j and m + z7j — > n, respectively p2[ . 



There is one final comment we wish to make. The precise identifica- 
tion that we have made between the anti-hermitian part of the Hamiltonian 
and the total rate for decays and inverse decays is intuitively appealing and 
useful for practical purposes. While ordinarily that allows us to determine 
the anti-hermitian part of the Hamiltonian by computing the relevant rates, 
there is one situation in which it does not. If the background contains neu- 
trinos, the rates for processes involving the background neutrinos depends 
on their number densities which, as already remarked, are modulated by the 
oscillation mechanism. It is then not obvious how to calculate the rates for 
such processes by the traditional way of calculating transition probability 
amplitudes. However, as a byproduct of our formalism we obtain the fol- 
lowing very precise procedure. We calculate the 21 and 12 elements of the 
self-energy matrix using the FTFT Feynman rules and then from Eqs. ( |119| ) 
and ( |121| ) obtain Td j. The dependence on the neutrino number densities 
shows up in the diagrams in which background neutrinos propagate in the 
internal lines of the loops. Thus, by using in those diagrams the same effec- 
tive neutrino propagator that was introduced in Section |^ to calculate the 
effects of the neutrino background on the dispersive terms, the effect of the 
modulation of the number densities by the oscillating amplitudes is properly 
taken into account in this case also. 
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Figure Captions 
Fig. 1. Neutrino background contribution to the self-energy. 
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